Local diffusion theory of localized waves in open media 

Chu-Shun Tian 1 , Sai-Kit Cheung 2 and Zhao-Qing Zhang 2 
1 Institut fur Theoretische Physik, Universitat zu Koln, D-50937 Koln, Germany 
Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong 

We report a first-principles study of static transport of localized waves in quasi-one-dimensional 
open media. We found that such transport, dominated by disorder-induced resonant transmissions, 
displays novel diffusive behavior. Our analytical predictions are entirely confirmed by numerical 
simulations. We showed that the prevailing self-consistent localization theory [van Tiggelen, et. 
al., Phys. Rev. Lett. 84, 4333 (2000)] is valid only if disorder-induced resonant transmissions are 
negligible. Our findings open a new direction in the study of Anderson localization in open media. 

PACS numbers: 42.25.Dd,71.23.An 
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Introduction. — In the past years experimental stud- 
ies of localization have been boosted due to the un- 
precedented level of manipulating ultracold atomic gases 
PQ, dielectric materials [2H7]> and elastic media [Ej- A 
key feature shared by many experimental setups [5H5] 
is that, there, one allows wave energies to leak out of 
systems through boundaries in order to facilitate mea- 
surements. Consequently, wave interferences interplay 
strongly with the wave energy leakage that, conceptually, 
enriches transport phenomena of localized waves while, 
technically, pushes forward developments of theoretical 
approaches. In particular, as a unique property of finite- 
sized samples, localized states in the sample center create 
resonant transmissions [10l[TT]. Although these transmis- 
sions are rare events, nevertheless, they contribute signif- 
icantly to average transmission. In fact, random matrix 
theory predicts that in quasi-one-dimensions (quasi-lD), 
the localization length measured from average transmis- 
sion can be four times larger than that measured from 
typical transmission [T2] . Recently, disorder- induced res- 
onant transmissions have found considerable practical 
applications. For example, they mimic a "resonator" 
with high-quality factors and thus are used to fabricate 
random laser [13] and to realize optical bistability [14] . 

However, to study this intriguing interplay has proved 
to be, in general, a formidable task. Confronting this 
challenge, a decade ago van Tiggelen, Lagendijk, and 
Wiersma took a bold stroke [TS], hypothesizing a so- 
called self-consistent local diffusion (SCLD) model for 
localization in open media (for a review, see Rcf. 16 ). 
They phenomenologically generalized the self-consistent 
localization theory of infinite media [16 by demanding 
the diffusion coefficient to be position-dependent so as 
to take boundary effects into account. The SCLD model 
(as well as its dynamic generalization [U [17]); having 
the advantages of physical transparency and method- 
ological simplicity over other approaches, is guiding con- 
siderable experimental and theoretical activities (e.g., 
Refs. [71 E □!)■ Yet, the validity of the SCLD model 
inside localized samples is largely unknown and, in fact, 
has been severely questioned by recent pulsed microwave 
experiments [7| . There, it was shown that the (dynamic) 



SCLD model [T7] fails to describe transport in quasi-lD 
localized samples at long times, when energies are mainly 
stored in long-lived modes. The dramatic discrepancy 
[7] between experimental measurements and theoretical 
predictions is conveying an opinion. That is, the highly 
non-local object of disorder-induced resonant transmis- 
sion |10j . which plays a decisive role in transport of lo- 
calized waves [llj . may not be captured by such a model. 

Motivated by these activities, we performed a first- 
principles study of static wave transport in quasi-lD lo- 
calized samples, i.e., L>^ with L and £ the sample and 
localization length, respectively. We predicted analyti- 
cally and confirmed numerically that in these systems, 
localized waves display a novel diffusion phenomenon. 
Our theory shows that the SCLD model is valid only 
if disorder-induced resonant transmissions are negligible. 
Our findings, capable of being generalized to higher di- 
mensions, may open a new direction in the study of An- 
derson localization in open media. 

Main results and qualitative discussions. — We consid- 
ered the wave intensity and particularly its spatial cor- 
relation function, y(x,x'). Our first-principles ana- 
lytic theory, justifying the static local diffusion equation: 
—d x D(x)d x y(x,x') = 8{x — x') with x (x') the distance 
of the observation (source) point from given boundary, 
leads to the following central results, (i) The local (or 
position-dependent) diffusion coefficient D(x) displays a 
novel scaling behavior. Specifically, D(x) depends on x 
via the scaling A = (L - x)x/{L£) [D = D(0)], 



D(x)/D =D oo (X). 

and the scaling function D QO (X) is ~ 



for A 



oo 



(1) 

(h) 



From (i) it follows that inside the sample, surprisingly, 
D(x) is enhanced drastically from the exponential decay, 

D(x)/D oce x2/(LS) e' x ^, £ < x < L/2 . (2) 

[D{x) = D(L - x) .} (iii) Eqs. and (J2J) are universal 
regardless of the time-reversal symmetry. Our results, 
while entirely confirmed by simulations (Fig. [I]), show 
that the SCLD model fails in localized samples. 
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FIG. 1: Comparing results obtained from numerical simu- 
lations, the analytical prediction S (solid lines) and the 
SCLD model (dashed lines). D{x)/Dq was computed numer- 
ically for two wave frequencies, cj = 1.65c/o (square) and 
oj = 0.72c/a (circle), and for five different sample lengths, 
L/C = 2.5, 5, 10, 15 and 20. 



Let us first present qualitative explanations of the main 
results (i)-(iii). The scaling behavior, Eq. (nj), finds its 
origin in wave interferences. Indeed, as waves pene- 
trate into a time-reversal medium, they may counter- 
propagate along the same loop and interfere with each 
other — the well-known weak localization [19]. However, 
different from infinite media, in open media, wave en- 
ergies leak out through the boundaries and as such, 
the probability (in the frequency domain) of forming a 
loop is finite in the static limit, which is oc [L — x)x/L 
in quasi-lD. Then, a one-loop wave interference correc- 
tion to D results, which is z-dependent and of the or- 
der of A = (L — x)x/(L£) . Furthermore, because A 
monotonously decreases from the sample midpoint to 
boundaries, as waves propagate towards the sample cen- 
ter, the propagation paths tend to form more (say n) 
loops with a probability oc \{L — x)x/L] n , leading to a 
wave interference correction ~ A™ . Thus, wave interfer- 
ences everywhere render D(x) depending on x via A . 

That the scaling behavior of local diffusion is far be- 
yond the reach of the SCLD model is best appreci- 
ated by broken time-reversal systems (unitary symme- 
try). There, the above one-loop wave interference is ab- 
sent and, therefore, the SCLD model ceases to work. [In- 
deed, for A — > , the linear term in the A-expansion of 
-Doo(A) now disappears.] Instead, two paths may take 
the same route, propagate in the same direction but visit 
individual scatterers at different times. As such, they 
form loops and equally contribute wave interference cor- 
rections to D . The perturbative A-expansion is thereby 
justified, with the leading order correction ~ A 2 . 

Having explained Eq. ([I]), let us estimate the asymp- 
totic form of .D 00(A) at A — > 00 . To this end we enjoy 
the universality of -Doo(A) and set L — ¥ 00 (semi- infinite 
media). In this simple case, on physical grounds, we ex- 



pect D(x)/D = £>oo(A) ~ e x lt for x > £ . Because 
of A = x/£, we find D M (A) ~ e~ x for A > 1 . This 
asymptotic form then gives Eq. (pi) for finite-sized sam- 
ples. Importantly, the significant enhancement from the 
exponential decay may be related to the fluctuation of 
the inverse localization length 7 in finite-sized samples. 
Indeed, for L ^> £ , the distribution of 7 is Gaussian, with 
the average and variance being £ _1 and 2/(£Z) , respec- 
tively [THr^D]. Averaging e~ lx , we obtained 



dje 



-^(-f-r 1 ) 2 



D(x) 



(3) 



for x deep inside the sample. 

Failure of SCLD model. — For simplicity we focus on 
classical scalar waves, and begin with testing the validity 
of the SCLD model [TS]. We performed numerical sim- 
ulations of the spatially-resolved wave intensity across 
a randomly layered medium, which is embedded in an 
air background and excited by a plane wave of (angular) 
frequency ui . The layer thickness is a , and the rela- 
tive permittivity at each layer fluctuates independently, 
with a uniform distribution in the interval [1 — u, 1 + o~] . 
Here a measures the degree of randomness of the system 
and throughout this work we considered non-reflecting 
boundaries. We set a — 0.7 and considered two wave 
frequencies, uj = 1.65c/a and 0.72c/a, where c is the 
speed of light in the air. We used the standard transfer 
matrix method to calculate the transmission coefficient, 
Tp , and wave intensity distribution, Ip (x) , for each con- 
figuration /3 . For each u>, we calculated the ensemble- 
averaged current j = (Tp) and wave intensity distribu- 
tion I(x) = (Ip (x)) of 2, 000, 000 realizations of dielectric 
disorders for different sample lengths. Since the current 
across the sample is uniform due to the conservation law, 
we used the relation: j = —D(x)d x I(x) to compute D(x) 
by presuming the static local diffusion equation. 

In Fig. [I] the results of D(x)/D obtained by simula- 
tions and by numerically solving the SCLD model are 
presented. First of all, they both show that D(x) tends 
to decay exponentially from the boundary in the limit: 
L — » 00 (dotted line) . The decay length (the localization 
length) £ was found to be the transport mean free path, 
which is 21a (50a) for u = IMc/a (0.72c/a). We there- 
fore rescale x into #/£ and present results for five different 
sample lengths, L/£ = 2.5, 5, 10, 15 and 20. We see that 
for different frequencies, the simulation results (squares 
and circles) overlap, signaling the scaling behavior inde- 
pendent of the parameters of random media. It is ob- 
vious that, except near the boundaries, the results from 
simulation are significantly larger than those from the 
SCLD model (dashed lines). The deviation is prominent 
for large L/£ , where the results from the SCLD model 
converge to the sum of two truncated exponentials, de- 
caying from their respective boundaries. Does localized 
waves in open media display diffusive transport? Our 
analytical prediction below provides a definitive answer 



to this conceptually important question. In particular, 
the analytical result of D(x)/Dq , Eq. pi), is in excellent 
agreement with numerical simulations (solid lines). 

Exact microscopic formalism. — Referring to a separate 
publication for technical details, we turn now to outline 
the proof and the strategy is as follows. In the framework 
of the supersymmetric field theory of localization [2Tll22| . 
we calculated explicitly the correlation function y(x, x') , 
and found that it solves the local diffusion equation. In 
doing so, we managed to calculate the weak localization 
correction, 5D(x) . Then, we found the Gell-Mann-Low 
equation of the local diffusion coefficient, D(x) = Dq + 
SD(x) , which eventually leads to Eq. pi). 

In the present context the supersymmetric technique 
has many advantages over others. The key ingredient is 
the introduction of a "spin" Q (to be defined below) that 
encapsulates wave interferences by fluctuations of the 
"spin direction" : the larger the fluctuation, the stronger 
the localization effect. With the help of the Q-spin, the 
picture of localization in open media is analogous to that 
of the more familiar problem — finite classical ferromag- 
netic spin chains (but formal treatments are not) . In the 
latter system, the two end spins are fixed and parallel, 
and the spin direction fluctuates elsewhere with a small 
(large) fluctuation amplitude closed to the chain ends 
(midpoint). Translated to the Q-spin language, such in- 
homogeneous fluctuations reflect the spatial inhomogene- 
ity of wave interferences in open media — the very mech- 
anism of local diffusion. Most importantly, the feature 
that the two Q-spins at the boundaries are "parallel" 
takes all disorder-induced resonant transmissions into full 
account. Thus, the phenomenology of local diffusion is 
substantiated by a completely microscopic formalism, al- 
beit in an elegant manner. 

Formally, the "spin" Q is an 8 x 8 supermatrix de- 
fined on the advanced-retarded (ar), bosonic-fermonic 
(bf) and time-reversal (tr) sector. The "ar" sector ac- 
commodates different analytic structures of the advanced 
(retarded) Green function. The "bf" sector accommo- 
dates the supersymmetry: the diagonal (off-diagonal) 
matrix elements are commuting (anti-commuting) num- 
bers. The "tr" sector accommodates the time-reversal 
symmetry. The leakage enters through the boundary con- 
straint: Q(0) = Q(L) = af . Here, o$ = diag(l,-l), 
X — ar, bf , tr . Then, y(x, x') is exactly expressed as 



y(x,x') 



D[Q]e 



-f[Q], 



xQ(x)(l-0(l-a*> 3 M Q(*')].(4) 



The action, F[Q] = -K^a j£ dxstr (d x Q) 2 [y the den- 
sity of states per unit length and "str" the supertrace), is 
the energy cost of the Q-field fluctuations. It introduces 
a characteristic scale-the localization length £ ex ttvDq . 
Importantly, the mean field Q(x) = erf r , compat- 
ible with the boundary constraint, minimizes the ac- 



tion and describes a vanishing wave intensity background 
across the sample. Bearing this in mind, we introduced 
the parametrization: Q = (1 + iW)af(l + iW)~ x , 
where W(x) anti-commutes with <jf r and vanishes at 
x = 0, L , and performed the W-expansion. By keep- 
ing the leading order W-expansions, we obtained from 
Eq. Q the bare correlation function &^(x, x') , which 
solves —D d x ?V(x, x') = S(x-x') with the boundary con- 
dition: <3f(0, x 1 ) = &(L,x') = 0. Calculating the leading 
wave interference corrections to &(x, x') , we found 

y(x, x') = &{x, x') -fdy &(x, y)d v 5D{y)d v W{y, x') (5) 

with the leading order weak localization correction read 



8D{x) <&{x,x) 1, 



W(x,x) 



(6) 



where ?V(x,x) = (L — x)x/(D L) , and a = — 1 for the 
orthogonal symmetry while a = for the unitary sym- 
metry. Eq. (p)]) justifies that D(x) depends on x via A . 

For the unitary symmetry the first term of Eq. Mm 
vanishes, reflecting that the one- loop interference is in- 
feasible. Thus, the local diffusion and scaling behavior 
of D(x) are universal concepts, extrinsic to the time- 
reversal symmetry that is required by the SCLD model. 

Scaling theory of local diffusion and disorder-induced 
resonant transmission. — We now make an important ob- 
servation: Q(X) = vD{x)/{£,\) and A , formally, play the 
role of the "Thouless conductance" [TH] and the "system 
size" , respectively. Indeed, from Eq. ffib we found 



d\nG 
dlnX 



= p(G) = -i + c 1 g- 



c 2 g~ 



, Q » 1 , (7) 



with C\ < for a — —I and c\ — 0, c 2 < for a = . 
This is fully analogous to the usual one-parameter scaling 
theory of quasi-lD localization [19] where, in particular, 
the perturbative expansion of the /3-function also finds 
its origin in weak localization [TBI |2"2] . 

From Eq. ([6]), we further found that the weak localiza- 
tion corrections of open and infinite media, including the 
coefficients, are identical except that the returning prob- 
ability, ty(x, x) , replaces that of infinite media. This 
duality persists in all the higher-order weak localization 
corrections and, as such, the /3-function here is identical 
to that of the usual one-parameter scaling theory [19j . 
Identifying this duality, we followed Refs. [19j [20] to ex- 
trapolate Eq. ^\ into the regime of Q -C 1 , obtaining 



0(G) =\nG, £«1 



(8) 



The scaling theory of local diffusion namely Eqs. ffh 
and (Isb is far beyond the reach of earlier theoretical stud- 
ies [22123] and has far-reaching consequences. (It may be 
reproduced within the formalism of Ref. |24j.) In particu- 
lar, Eq. B gives G(X) ex e~ A for A ->• ex) and thus Eq. (B 
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FIG. 2: Simulations results obtained by eliminating the high- 
transmission (T > T c ) states in the original ensemble (uj = 
1.65c/a and L/£ = 10). The source is placed at x — . 



which, as shown below, fully captures the rare disorder- 
induced resonant transmission. Solving the static local 
diffusion equation, we found that the ensemble-averaged 
transmission is (T(L)) ex [/„ dx/D(x)}~ 1 (This is re- 
gardless of the explicit form of D(x) , as first noticed 
in Ref. [IS].) Inserting Eq. (2) into it gives (T(L)) oc 
e -V( 4 ?) f or l ^ £ On the other hand, noticing that 
T(L) — e~~ fL , we found that the typical transmission 
gives d(\nT(L))/dL = — £ _1 from the Gaussian distri- 
bution of 7 [cf. Eq. pb]. Thus, the localization lengths 
obtained by the arithmetic and geometric means differ by 
a factor 4 irrespective of orthogonal or unitary symme- 
try. This is in agreement with the result of the random 
matrix theory, and is because that (T(L)) is dominated 
by disorder-induced resonant transmissions |12j . 

To further study effects of rare high-transmission states 
we analyzed all the samples in Fig.[l|(w = 1.65c/a). First 
of all, the distribution of InT is well fitted by the nor- 
mal distribution, with the average w —L/£ and the vari- 
ance ss 2L/£ . This confirms the Gaussian distribution of 
7 . Then, we intentionally eliminated a small fraction of 
high-transmission (T > T c ) states (composed mostly of 
singly localized states and of a small portion of necklace 
states [25]), from the original ensemble (L/t; = 10) and 
re-computed D(x)/Do . As shown in Fig. ^1 strikingly, 
even when the fraction of removed states is as small as 
0.6% (solid line, in red), the simulation result deviates 
drastically from the original one (solid line, in black) and 
is asymmetric, signaling the breakdown of local diffusion. 
Thus, we found that rare high-transmission states are es- 
sential to establish local diffusion and scaling behavior. 

Conclusions. — We found in quasi- ID localized samples 
a scaling behavior of the (static) local diffusion coefficient 
capturing all the rare disorder-induced resonant trans- 



mission. Our findings show unambiguously that the pre- 
vailing SCLD model is valid only if rare disorder-induced 
resonant transmissions are negligible which, nevertheless, 
play a decisive role in transport of localized waves. 

The found phenomenon is intrinsic to finite-sized sam- 
ples with open boundaries and does not exist in an in- 
finite sample. It is an unconventional diffusion phe- 
nomenon in the sense that the diffusion coefficient can 
drop by many orders of magnitude as the position 
changes from the boundary to the midpoint. (For dif- 
fusive samples, such a position-dependence is weak and 
thus does not lead to any interesting phenomenon other 
than ordinary diffusion.) It is such a drastic change (in 
the diffusion coefficient) that leads to a global localization 
behavior as shown in the scaling of the average transmis- 
sion which decays exponentially with sample size. 

Our theory has many immediate applications. For ex- 
ample, it can be directly used to study the speckle pat- 
tern of scattered waves which has recently attracted con- 
siderable attentions. It may also be generalized to higher 
dimensions for studying disorder-induced resonant trans- 
missions in (or close to) the localization regime. This 
issue has practical applications such as random laser. 

We thank A. Z. Genack for discussions, and A. Altland, 
D. Basko, A. A. Chabanov, M. Garst, A. Kamenev, and 
T. Nattermann for conversations. Work supported by 
SFB TR12 of the DFG and the HK RGC No. 604506. 



[1] 



J. Billy et ai, Nature 453, 891 (2008); G. Roati et ai, 

ibid. 453, 895 (2008). 
[2] T. Schwartz, et ai, Nature 446, 52 (2007); 
[3] K. Yu. Bliokh et al., Phys. Rev. Lett. 97, 243904 (2006). 

C. Toninelli, et ai, Phys. Rev. Lett. 101, 123901 (2008). 

A. A. Chabanov, et. ai, Nature 404, 850 (2000). 

M. Stozer, et ai, Phys. Rev. Lett. 96, 063904 (2006). 

Z. Q. Zhang et ai, Phys. Rev. B 79, 144203 (2009). 

H. Hu, et al, Nature Physics 4, 945 (2008). 

A. A. Chabanov, Z. Q. Zhang, and A. Z. Genack, Phys. 
Rev. Lett. 90, 203903 (2003). 

[10] M. Y. Azbel, Phys. Rev. B 28, 4106 (1983). 

[11] K. Y. Bliokh, et. ai, Rev. Mod. Phys. 80, 1201 (2008). 

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

H. Cao, et. ai, Phys. Rev. Lett. 82, 2278 (1999). 

I. V. Shadrivov, et. al, Phys. Rev. Lett. 104, 123902 

(2010). 

B. A. van Tiggelen, A. Lagendijk, and D. S. Wiersma, 
Phys. Rev. Lett. 84, 4333 (2000). 

P. Wolfle and D. Vollhardt, to appear in Fifty Years of 
Anderson localization, ed. E. Abrahams (World Scien- 
tific, Singapore, 2010). 

S. E. Skipetrov and B. A. van Tiggelen, Phys. Rev. Lett. 
92, 113901 (2004); 96, 043902 (2006). 

[18] N. Cherroret, et al, Phys. Rev. B 80, 037101 (2009). 
E. Abrahams, et al, Phys. Rev. Lett. 42, 673 (1979). 
P. W. Anderson, et ai, Phys. Rev. B 22, 3519 (1980). 

C. Tian, Phys. Rev. B 77, 064205 (2008). 
K. B. Efetov, Adv. Phys. 32, 53 (1983). 



[4] 
[5] 
[6] 
[7] 
[8] 
[9] 



[12] 
[13] 
[14] 

[15] 

[16] 

[17] 



[19] 
[20] 
[21] 

[22] 



[23] N. Cherroret and S. E. Skipetrov, Phys. Rev. E 77, and Perspectives eds. W. Janke and A. Pelster (World 

046608 (2008). Scientific, Singapore, 2008). 

[24] P. W. Brouwer and A. Altland, Phys. Rev. B 78, 075304 [25] J. B. Pendry, Adv. Phys. 43, 461 (1994). 
(2008); A. Altland et. al, in Path Integrals-New Trends 



